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A formalism for studying spontaneous decay of an excited 
two-level atom in the presence of dispersing and absorbing 
dielectric bodies is developed. An integral equation, which 
is suitable for numerical solution, is derived for the atomic 
upper-state-probability amplitude. The emission pattern and 
the power spectrum of the emitted light are expressed in terms 
of the Green tensor of the dielectric-matter formation includ- 
ing absorption and dispersion. The theory is applied to the 
spontaneous decay of an excited atom at the center of a three- 
layered spherical cavity, with the cavity wall being modeled 
by a band-gap dielectric of Lorentz type. Both weak coupling 
and strong coupling are studied, the latter with special em- 
phasis on the cases where the atomic transition is (i) in the 
normal-dispersion zone near the medium resonance and (ii) 
in the anomalous-dispersion zone associated with the band 
gap. In a single-resonance approximation, conditions of the 
appearance of Rabi oscillations and closed solutions to the 
evolution of the atomic state population are derived, which 
are in good agreement with the exact numerical results. 

PACS numbers: 42.50.Ct, 12.20.-m, 42. 60. Da, 42.50.Lc 



I. INTRODUCTION 

It is well known that the spontaneous decay of an ex- 
cited atom can be strongly modified when it is placed 
inside a microcavity . There are typically two quali- 
tatively different regimes: the weak-coupling regime and 
the strong-coupling regime. The weak-coupling regime 
is characterized by monotonous exponential decay, the 
decay rate being enhanced or reduced compared to the 
free-space value depending on whether the atomic tran- 
sition frequency fits a cavity resonance or not. The 
strong-coupling regime, in contrast, is characterized by 
reversible Rabi oscillations where the energy of the ini- 
tially atom is periodically exchanged between the atom 
and the field. This usually requires that the emission 
is in resonance with a high-quality cavity mode. Re- 
cent progress in constructing certain types of microcavi- 
ties such as microspheres has rendered it possible to ap- 
proach the ultimate quality level determined by intrinsic 
material losses || , so that the question of the influence 
of absorbing material on spontaneous decay has been of 
increasing interest. 

Effects of material losses on the lifetime of an excited 
atom have been studied within Fermi's golden-rule ap- 
proach |i|-[TT|. In Q the mode structure of a micro- 
sphere without and with absorber dopant atoms, which 



is modeled by a constant and a Lorentzian dielectric func- 
tion respectively, is considered. The spontaneous emis- 
sion rate and the radiation intensity as a function of the 
atomic transition frequency is examined in [p"3[ for an 
atom in a Fabry-Perot cavity filled with a Lorentz-type 
dielectric in the case of strong medium-cavity interaction 
but weak atom-field interaction. 

In this paper we present a theory of the spontaneous 
decay of an excited two-level atom in the presence of 
arbitrary dispersing and absorbing dielectric bodies. We 
apply the theory to the spontaneous decay of an atom in a 
spherical microcavity of given complex-valued refractive- 
index profile, as it is typically the case in experimental 
implementations. The formalism enables us to include in 
the theory absorption and dispersion in a consistent way 
and to give a unified treatment of spontaneous emission, 
without restriction to a particular coupling regime. 

The plan of the paper is as follows. In Section H, a 
recently developed quantization scheme for the electro- 
magnetic field in the presence of dispersing and absorbing 
dielectric bodies [pd|-p"7| is extended in order to include 
in the theory the resonant interaction of the field with 
a two-level atom. From the Hamiltonian of the com- 
posed system, an integral equation governing the tem- 
poral evolution of the upper-level-probability amplitude 
of the atom is derived, the integral kernel being deter- 
mined by the Green tensor of the classical, phenomeno- 
logical Maxwell equations for the dielectric-assisted elec- 
tromagnetic field. General expressions for the emission 
pattern and the power spectrum are derived in terms of 
the atomic parameters and, via the Green tensor, the 
cavity parameters of the dielectric-matter configuration. 

In Section III, the theory is used to examine the spon- 
taneous decay of an excited two-level atom inside a spher- 
ical cavity, with special emphasis on the intrinsic disper- 
sion and absorption of the wall material. Spherical mi- 
crocavities have been very attractive systems for both 
fundamental research in cavity quantum electrodynam- 
ics and applications in optoelectronics (see, e.g., |5 18 - 2C|] 
and references therein). Changes in the level shifts and 
lifetimes of atoms inside or near the surfaces of spherical 
micro-structures have been studied theoretically jglj , the 
latter also experimentally |22| , and results on the strong- 
coupling regime have been reported [p3|-p5| . The theoret- 
ical results have been typically based on standard mode 
expansion, which fails for intrinsically dispersing and ab- 
sorbing material. Here, we consider a spherical three- 
layered structure, where the middle layer is assumed to be 
a (single-resonance) band-gap dielectric of Lorentz type, 
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whereas the outer and inner layers are vacuum. The 
strength of the atom-field coupling, which is essentially 
determined by the imaginary part of the Green tensor at 
the position of the atom, is analyzed and the positions, 
heights, and widths of the possible cavity resonances 
are calculated. In particular, the most pronounced res- 
onances are observed within the band gap, the widths 
of which are proportional to the intrinsic damping con- 
stant of the medium. In a single-resonance approxima- 
tion, conditions of the strong-coupling regime and closed 
expressions for the atomic upper-state-population proba- 
bility are derived, which are in good agreement with the 
numerical solutions. Finally, a summary and conclusions 
are given in Section IV. 



II. GENERAL FORMALISM 

A. Quantization scheme 

Let us first consider the electromagnetic field in the 
presence of dispersing and absorbing dielectric bodies 
without additional atomic sources. Following Jl^-IT^I, we 
represent the electric-field operator E in the form 

E(r) =E(+)(r)+E(-)(r), EJ(->(r) = [E (+) 00] f (1) 

E' + '(r) = / dwE(r,w), (2) 
Jo 

and the magnetic-field operator B accordingly. The op- 
erators E and B then satisfy the Maxwell equations 



V ■ B(r,w) = 0, 
V- e e(r, u)E(r, oj) = p(r,u), 
V x E(r,w) = wjB(r,u;), 



(3) 
(4) 
(5) 



V x B(r, u) = -i 3 e(r, u/)E(r, u) + /x j(r, w), (6) 

where the complex permittivity e(r, to) is a function of 
frequency and space, the real part (er) and the imagi- 
nary part (ei) of which satisfy (for any r) the Kramers- 
Kronig relations. The operator noise charge and current 
densities p(r, lS) and j(r,w) respectively, which are asso- 
ciated with absorption, are related to the operator noise 
polarization P(r, ui) as 



p(r,u) = -V-P(r,w), 
j(r,w) = -iuP(r,u), 



where 



with 



E(r,w) 



—ei(r,w) f(r, lu) 
ir 



(7) 
(8) 

(9) 



Mr,u;)JJ(r\Lu') = %«5(r-r')%-u/), (10) 



/i(r,w),/ J (r' ) w / ) 



= = 



f}(r,u),f (r>,u>) . (11) 



From Eqs. (Q) - (||) it follows that E can be written in 
the form 



E(r,«)=*4/ 7 



dV Vei(r',^)G(r, r', to) ■ f (r', w), (12) 



and B = (iui) 1 V x E accordingly, where G(r,r',ui) is 
the classical Green tensor satisfying the equation 



r, ,2 



■ e(r, w) — V x Vx 



G(r,r',u;) = -8{r-r') (13) 



together with the boundary condition at infinity [S(r) 
is the dyadic <5-function] . In this way, the electric- and 
magnetic-field strengths are expressed in terms of a con- 
tinuum set of bosonic fields f and F, which play the 
role of the fundamental (dynamical) variables of the com- 
posed system (electromagnetic field and the medium in- 
cluding the dissipative system) whose Hamiltonian is 



H = J d 3 r dwhu}&(T,w) ■ f(r,w). 



(14) 



Using Eqs. ( |l2| ) [together with Eqs. ([[]) and (g)], we can 
also express the scalar potential (p and the vector poten- 
tial A of the electromagnetic field in terms of the fun- 
damental bosonic fields. In particular, in the Coulomb 
gauge we obtain 



A(r) 



-V0(r)=E"(r), 
duj (iw) _1 E (r,w) + H.c. 



where 



E^(r) = J d 3 r' S x ^(r - r') • E(r'), 



(15) 
(16) 

(17) 



with S ± (r) and S"(r) being the transverse and longitu- 
dinal 5-functions respectively. 

We now consider the interaction of the medium-assis- 
ted electromagnetic field with additional point charges 
q a . Applying the minimal-coupling scheme, we may write 
the complete Hamiltonian in the form 



H = 



J d 3 v dufkJ^(r,u})-f(r,Lj) 



E 



2m r 



p a - q a A(r a ) ■ p a - q a A(r a ) 



+ ^ I d 3 rp A (r)£ A (r) + / d 3 r p A (r)^(r), (18) 
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where f Q is the position operator and p a is the canonical 
momentum operator of the ath charged particle of mass 
m a . The Hamiltonian (jlq) consists of four terms. The 
first term is the energy (JlJ) observed when the particles 
are absent. The second term is the kinetic energy of the 
particles, and the third term is their Coulomb energy, 
where the potential (pA can be given by 



with 



<Mr) 



PA ( 



dr' ■ 



PA(r') 



47T£o r — r' 



q a S(r - r a ) 



(19) 



(20) 



being the charge density. The last term is the Coulomb 
energy of interaction of the particles with the medium. 
Note that all terms are expressed in terms of the dynam- 
ical variables f (r, w), f'(r, w), r a , p a . 



B. Dynamics of an excited two- level atom 

Let us consider a neutral two- level atom (position r A , 
transition frequency uja) that resonantly interacts with 
radiation via an electric-dipole transition (dipole moment 
fi). In this case, the electric-dipole approximation and 
the rotating wave approximation apply, and the minimal- 
coupling Hamiltonian ( |l8| ) simplifies to (Appendix A) 

H = J d 3 r duTiujP(r,u)-f(r,uj) 

+ ±huj A a z - [«r t E(+) (r A ) ■ fi + H.c.l , (21) 



where a, a\ and a z are the Pauli operators of the two- 
level atom. 

When the atom is initially in the upper state and the 
rest of the system is in the vacuum, then the system wave 
function at time t can be written as 

\m) = c u {t)e-^i^\ U )m) 

poo 

d 3 r J ^C /4 (r, W ,i)e- l( ^ A/2)t |0|{l,(r,^}), (22) 

where \u) and |^) respectively are the upper and lower 
atomic states, |{0}) is the vacuum state of the rest of the 
system, and |{li(r,o>)}) is the state, where the latter is 
excited in a single-quantum Fock state. Here and in the 
following we adopt the convention of summation over re- 
peated vector-component indices. The Schrodinger equa- 
tion yields 



C u {t) = - 



dui 



d A v 



\Ai(r,w) 



x G^(r A , r, u) C H (r, u, i)^"-"*)* , (23) 



Cu(r,u,t) = -^=^ V £ i( r ^) 

xG* i (r A ,r,o;)C„(t) e i (™)*. 



(24) 



We now substitute the result of formal integration of 
Eq. (|J) [C H {r,Lu,0) = 0] into Eq. (f|). Making use of 
the relationship 



Im G fci (r, t',uj) 



f LlI 2 

/ d3s ^2" ei(s,u)Gk m (r,s,uj)Gl m (r' ,s,uj) 
we obtain the integro-differential equation 



C u (t) = / dt 1 K{t-t')C u (t'), 
Jo 

with the kernel function 



(25) 



(26) 



K(t - t') = 



x / duje-^-^^'hmGi^r^Tx^) (27) 



(fc A = loa/c). In the spirit of the rotating wave approxi- 
mation used we have set uj = w A in the integral in Eq. (|27]) . 

Taking the time integral of both sides of Eq. (p6|), we 
easily derive, on changing the order of integrations on the 
right-hand side, 



C u (t) = / dt'K(t 
Jo 



t')C u {t') + 1 



(28) 



[C„(0) -1], where 
fc A /ij/ij 



K(t - t 1 ) 



< / duj 
la 



ImG,,(rA,rA,a)) 
i(oj - uj a ) 



-i(u — uJA.)(t—t ) 



- 1 



(29) 



The integral equation (|2q ) is a well-known Volterra equa- 
tion of the second kind. An algorithm for solving such an 
integral equation numerically can be found, e.g., in |p6| . 
It is worth noting that the integro-differential equation 
( p6| ) and the equivalent integral equation ( |28] ) apply to 
the spontaneous decay of an atom in the presence of an 
arbitrary configuration of dispersing and absorbing di- 
electric bodies. All the matter parameters that are rele- 
vant for the atomic evolution are contained, via the Green 
tensor, in the kernel functions (|27j ) and ((2^). In particu- 
lar when absorption is disregarded and the permittivity 
is regarded as being a real frequency-independent quan- 
tity (which of course can change with space), then the 
formalism yields the results of standard mode decompo- 
sition (see, e.g. @-§9|). 

When the Markov approximation applies, i.e., when in 
a coarse-grained description of the atomic motion mem- 
ory effects are disregarded, then we may let 
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J(u A -u)(t-t') _ 

i{u> A - oj) 



(30) 



inEq. @ l((x) -- 
value], and thus 



:ttS(x) + iP/x; V denotes the principal 



where 



K(t-t') = -UA-iSto). 



A = ^%^IrnG y (r A ,r A , WA ) 
tier, 



and 



s ^ _ 2k\^j f 00 ^ Im G t j(r A ,r A ,ui) 



u> — uj a 



(31) 



(32) 



(33) 



Substituting into Eq. ( p8| ) for the kernel function the ex- 
pression (|3l|), we obtain the familiar result that 



C u (t) = exp[-±(A-i6w)t] 



(34) 



Obviously, this result is also obtained if in the integral in 
Eq. (|2^) C u (t') is replaced by C u (t) and then the integral 
is approximated by Q{loa —uj). Note that the expressions 
( |32| ) and (^) for the decay rate and the line shift, re- 
spectively, are in full agreement with the results in |Tc| ] . 

It is well known that the Markov approximation is an 
excellent approximation for describing the radiative de- 
cay of an excited atom in free space. In order to study 
the case where the atom is surrounded by dielectric mat- 
ter, we assume that the atom is localized in a more or 
less small free-space region, so that the Green tensor at 
the position of the atom reads ]3(| 

G(r A ,r A ,w) = G v (r A ,r A ,cj) + G R (r A , r A , w), (35) 

where G v is the vacuum Green tensor, with 



v ^ 
ImG (r A ,r A ,w) = - — I 
one 



(36) 



(for the vacuum Green tensor, see, e.g., |||3Tj), and G R 
describes the effects of reflections at the (surfaces of dis- 
continuity of the) surrounding medium. The contribu- 
tion of G v to K can then be treated in the Markov ap- 
proximation. Application of Eqs. (31) - ( |33| ) yields the 
well-known vacuum decay rate 



A 



J? 3 



AM 



3?i7reo 



(37) 



and a divergent contribution to the vacuum Lamb shift 
which may be omitted, since the (renormalized) vacuum 
Lamb shift may be thought of as being included in the 
atomic transition frequency uj a . In this way, Eq. ( p9| ) 
takes the form 



K{t-t') =-±Ao + 



k A ^i^j 

hirer, 



Im Gg(r A ,r A , uj) 
- uj a ) 



(38) 



The integral equation ( p8| ) together with the kernel func- 
tion (|8|) can be regarded as the basic equation for study- 
ing the influence of an arbitrary configuration of dispers- 
ing and absorbing dielectric matter on the spontaneous 
decay of an excited atom. 



C. Emission pattern 

The intensity of the light registered by a point-like pho- 
todetector at position r and time t is given by 



7(r,i) = (V^)|E(-)(r).E(+)(r)|V>(*)> 



(39) 



To obtain the emission pattern associated with the spon- 
taneous decay of an excited atom in the presence of 
dispersing and absorbing dielectric matter, we combine 
Eqs. (0), (§), @, and @. After some algebra we de- 
rive, on using Eq. (Ea), 



TTCo Jo 



d£ 



C u (t') 



< I dujlmGij(r,r A ,uj)e 
Jo 



— i(cd— oja)(*— t ) 



(40) 



where we have again set u> — oj a in the frequency integral. 
Again, all the relevant dielectric-matter parameters are 
contained in the Green tensor. In contrast to Eq. ( |2S| ) 
together with the kernel function (^8|), Eq. ([foj) requires 
information about the Green tensor at different space 
points. In particular, its dependence on space and fre- 
quency essentially determines the retardation effects. 

In the simplest case of the atom being in free space we 
have 



Hjlm Gjj (r, r A 



1 



8iirp 



x ( e iup / c — e~ iujp / c 



PP P 

P 2 



(41) 



(p = r — r A ). We substitute the expressions ( |34| ) (with A 
— Aq) and ( fi"l|) into Eq. (|40|), calculate the time integral, 
and extend in the frequency integral the lower limit to 
— oo, which then can be calculated by contour integra- 
tion, 



duo (e wp l c - e ~ Mp/c ) 



3-(A /2+iuj' A )t _ p —iuit 



i[u - K - M /2)] 
= -27rcxp[(-iA - icj' A ){t - p/c)] 0(t - p/c), (42) 
where 

uj' a = uj a — TjSui (43) 

[©(a;), unit step function]. We thus derive the well-known 
(far-field) result that 



I(r,f) 



/ fc AM smfl \ 2 e -A 0(t - P /c) Q(t _ /c) (44) 
V Aireop J 
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where 9 is the angle between p and p. 

Let us return to the general expression (fhij). If retar- 
dation is ignored and the Markov approximation applies, 
then we can replace, for all r, C u (t') by C u (t) in the time 
integral in Eq. ( f4(i| ) and approximate the time integral by 
C(uja — w). We obtain 



J(r,i) = |F(r,r A ,w A )| 2 e- 



-At 



(45) 



where 



Fj(r, t a ,wa) 



7re 



/ dwImGy(r, rA,w)C(wA - w). 
Jo 



(46) 




FIG. 1. Scheme of the spherical cavity. 



D. Emitted-light spectrum 

Next, let us consider the time-dependent power spec- 
trum of the emitted light, which for sufficiently small 
passband width of the spectral apparatus can be given 
by (see, e.g., ||) 



S(t,ub,T) = I dt 2 [ dh 



-iui a (t 2 -t x ) 



JO 

x (E<->(r,ta). £(+)(!■,*!)) 



(47) 



where u>s is the setting frequency of the spectral appara- 
tus and T is the operating-time interval of the detector. 
In close analogyto the derivation of Eq. (^) , we combine 
Eqs. (0), (|), (£!), and (§|) and use the relation (§|) to 
obtain 



S(r, Ws ,T) = ]T 



Treo Jo 



i(UB-UA)h dt'C u (t') 



duj Im Gij (r, ta , u)e 



-i(w-w A )(tl-t') 



(48) 



Further calculation again requires knowledge of the 
Green tensor of the problem. 

Let us apply Eq. (|4g) to the free-space case. Following 
the line that has led from Eq. ([h]) to Eq. Q), we find 
that 



S(r,ws,T) = 



47T60P 



,[-Ao/2+%(us-<»a)](T-p/c) _ i 



lu s - lu' a + iA /2 



®{T-p/c). (49) 



In particular for T — * oo, we recognize the well-known 
Lorentzian: 



lim S(r,u; s ,T) = 

T — >oo 



k\ fi sin 9 
1 



[c^s - (c^a - 5u/2)f + All 'A 



(50) 



When retardation is ignored and the Markov approxi- 
mation applies, then Eq. ( [48| ) can be simplified in a sim- 
ilar way as Eq. (|40|). In close analogy to the derivation 
of Eq. (Ha) we may write 



S(r,uj s ,T) = \F(r,r A ,o JA )\ 2 

e {-A/2+i[uj a -(u A -&u/2)}\T 



1 



lu s - (cj a - 5ai/2) + iA/2 
with F(r,rA,^A) being defined by Eq. (M) 



(51) 



III. APPLICATION TO A SPHERICAL CAVITY 
A. The model 

We apply the formalism developed in Sec. || to the 
spontaneous decay of an excited two-level atom placed 
inside a spherical three-layered structure (Fig. [l]). The 
outer layer (r > i?i) and the inner layer (0 < r < R2) 
are vacuum, while the middle layer (R2 <r<R\) is a 
dispersing and absorbing dielectric. The Green tensor of 
the configuration is given in Appendix [B]. 

We have performed the calculations assuming a 
Lorentz-type dielectric with a single resonance (in the 
relevant frequency region): 



e(w) = 1 + 



(52) 



Here, up is the plasma frequency, which is proportional 
to the square root of the number density of the Lorentz 
oscillators and plays the role of the coupling constant be- 
tween the medium polarization and the electromagnetic 
field, and lot and 7, respectively, are the position and 
the width of the medium resonance. The plots in Fig. ^ 
of the real and imaginary parts of the index of refraction 



i(u) = y t{io) = ur(uj) + ini{uj) 



(53) 



illustrate a typical band-gap behavior of the configura- 
tion. 
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FIG. 2. (a) The real part 71r(u;) and (b) the imaginary 
part n\(u}) of the complex refractive index are shown for lop = 
0.5 lot and 7= 10 -4 ojt (solid line), j = 10~ s lot (dashed line), 
and 7 = 10~ 2 lot (dotted line). The longitudinal frequency lol 



y/u> T +LOp is Wl ! 



: 1.12 Wxi an d hence u>l — wt — 0.12 wt- 



B. Weak-coupling regime 

In the weak-coupling regime, the excited atomic state 
decays exponentially [Eq. (p4|)]. For simplicity let us as- 
sume that the atom is positioned at the center of the 
cavity. From Eqs. @, @, @), and ( |B22| ), the cavity- 
modified decay rate is then found to be 



A = A{lo a )A q , 
where Aq is decay rate in free space, Eq. (|3 
A(u) = l + RcCfH , 



and 



(54) 



(55) 



with C% 3 (w) being given by Eqs. (|b|) - flB19p . Note that 
if mode expansion applies, A(lo) would correspond to the 
change of the density of modes due to the presence of the 
cavity. 

When far from the medium resonance absorption is 
disregarded and hence the (frequency-independent) re- 
fractive index is assumed to be real, then previous re- 
sults obtained by mode decomposition can be recovered. 
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FIG. 3. The function A(uo), Eq. (B5|), is shown for R 2 = 
30 At, c?=At, up=0.5wt, and 7=10~ lot- The curves in the 

inset correspond to 7 = 10~ 2 lot (solid line), 7 = 2 x 10 -2 lot 

-2 



(dashed line), and 7 



= 10- 

5 x 10" 



' lot (dotted line) . 



In particular, for R2 — > we recognize the decay rate 
obtained in [^l],|5) for microspheres and liquid droplets. 
However, it should be pointed out that even far from the 
medium resonance the imaginary part of the refractive 
index cannot be set equal to zero in general, since the 
contribution to the decay rate of the nonradiative decay 
associated with absorption increases ~ R^ 3 for decreas- 
ing i?2 (and nonvanishing imaginary part of the refractive 
index) @. 

Let us restrict our attention to a true microcavity 
(R2LVA/ 'c^> 1). From Fig. || it is seen that the rate of 
spontaneous decay sensitively depends on the transition 
frequency. Narrow-band enhancement of spontaneous de- 
cay (A > 1) alternates with broadband inhibition (A< 1). 
The frequencies where the maxima of enhancement are 
observed correspond to the resonance frequencies of the 
cavity. Within the band gap the heights and widths of the 
frequency intervals in which spontaneous decay is feasible 
are essentially determined by the material losses. Out- 
side the band-gap zone the change of the decay rate is less 
pronounced because of the relatively large input -output 
coupling, the (small) material losses being of secondary 
importance. 

When R 2 u j/c^ > 1 and exp[— n\{R\ — R2)u)/d\ -C 1, then 
Eqs. |3|) - ( p!9| ) drastically simplify and A(lj) [Eq. " 
reads 



A(u>) ~ Re 



n(uj) — i tan(i?2<x>/c) 
1 — in(u>) tan(i?2w/c) 



(56) 



The positions Lo m of the maxima of A(ui) can then be 
obtained from the equation dA/dtu = 0. As long as n{io) 
can be regarded as being slowly varying compared with 
tan(a;i?2/c), we may neglect dri/dut and thus determine 
the resonance frequencies from the equation 

2ni(u) m ) tan(i? 2 w m /c) 

~ |n(w m )| 2 -l-^(|n(o; m )| 2 -l) 2 + 4n2( Wm ). (57) 



G 



Note that Eq. (|57]) is exact when it is regarded as condi- 
tional equation of R2 for a desired resonance frequency. 

In the band-gap zone we may assume that m 3> riR (see 
Fig. ||) . From Eqs. (|56|) and ( p7| ) it then follows that the 
maximum values A(to m ) and half widths at half maxi- 
mum, 6u> m , of the cavity resonance lines, i.e., the regions 
where enhanced spontaneous decay can be observed, are 
given by 



A(u> Tl 



i\{iO m ) + 1 2 V / K-^)«-^) 



«R(w m ) 



7W„ 



duJr, 



R2A(uj m ) 



(58) 



(59) 



where we have assumed that the material losses are small, 
i.e., 7<Cwt, wp, Wp/wx- Equations (p8|) and (59) reveal 



that in the approximation made the heights (widths) of 
the resonance lines (Lorentzians) are proportional (in- 
versely proportional) to 7, the highest and narrowest line 
being in the center of the band gap. Note that the prod- 
uct A(u> m )6uj m does not depend on 7. 

Outside the band-gap zone the inequality tt-r ^> n\ is 
typically valid (see Fig. |^), and thus Eqs. (^) and (|57 
yield 



A(w m ) ~ n R (w 



R2A{uj m ) 



(60) 



(61) 



The heights and widths of the resonance lines are now 
seen to be (approximately) independent of 7. The lines 
become higher and narrower if uo m becomes close to ljt . 

The widths of the resonance lines are responsible for 
the damping of intracavity fields. There are two damping 
mechanisms: photon leakage to the outside of the cavity 
and photon absorption by the cavity-wall material. From 
the analysis given above it is seen that the first mecha- 
nism is the dominant one outside the band gap where 
normal dispersion (dna/dw > 0) is observed, while the 
latter dominates inside the band gap where anomalous 
dispersion {dn-^/duj < 0) is observed. To illustrate this in 
more detail, we have calculated the amount of radiation 
energy observed outside the cavity, 



W = 2ce 



dt 



/>7T 

> / dO p 2 sin 9 J(r, t) 
Jo 



(62) 



(p > i?i), and compared it with the emitted energy in 
free space Wo = Hu>a- Assuming without loss of gener- 
ality that the atomic transition dipole is z-oriented and 
restricting our attention to the relevant far-field contri- 
bution, from Eqs. (|45|) and ( |46| ) together with Eq. (Bl) 
and Eqs. ( B2S ) - ( jB25| ) we derive (see Appendix y) 




C0a/(Ot 

FIG. 4. The amount of radiation energy W, Eq. (|63[), 
observed outside the cavity is shown as a function of the 
atomic transition frequency lja for 7 = 1CP 2 ljt (solid line), 
7 = 2 x 1CT 2 lot (dashed line), and 7 = 5 x 1CT 2 lot (dotted 
line). The other parameters are the same as in Fig. H. 



W \A%^1 
Wo l + HeCf(o;A) 



(63) 



Examples of the dependence on the atomic transition fre- 
quency of the amount of radiation energy observed out- 
side the cavity are plotted in Fig. ||. It is seen that inside 
the gap most of the energy emitted by the atom is ab- 
sorbed by the cavity wall in the course of time, while 
outside the gap the absorption is (for the chosen values 
of 7) much less significant. Note that with increasing 
value of 7 the band gap is smoothed a little bit, and thus 
the fraction of light that escapes from the cavity can in- 
crease. 



C. Strong-coupling regime 

The strength of the atom-field coupling increases when 
the atomic transition frequency loa approaches a cavity- 
resonance frequency u) m . In order to gain insight into the 
strong-coupling regime, let us first consider the limiting 
case of one cavity-resonance line being involved in the 
atom-field interaction. Indeed, when wa is close to u> m , 
then contributions from the other resonance lines become 
small. Using Eqs. (|32|), (|5^), and ([slf), and recalling that 
A(to) behaves like a Lorentzian in the vicinity of u> m , we 
may simplify Eq. (|2^) to 

K(t - t') ~ -p- AMiSiOmfe- 1 ^-^-^ 



2tt 



duj 



i(uj—u m )(t — t') 



-\A A{ujm)5Lo m e 



(lo - LU m ) 2 + (Suj m ) 2 



(64) 
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FIG. 5. The temporal evolution of the occupation prob- 
ability \C u (t)\ 2 of the upper atomic state is shown for R2 = 
30 A T , d = At, wp=0.5w T , w A = 1.046448w T , v4 A T /(2c) = 
10 -6 , and 7 = 10~ 4 ojt (solid line), 7 = 5 x 10~ 4 ^t (dashed 
line), 7 = 10~ ut (dotted line). For comparison, the expo- 
nential decay in free space is shown (dashed-dotted line). 



Substituting this expression into Eq. (26) and differenti- 
ating both sides of the resulting equation with regard to 
time, we arrive at 

C u (t) + [<(w m -w A ) + Su m ] C u (t) + (fl/2) 2 C u (t) = 

(65) 



where 



y 1 2A A(u] m )5u> r , 



(66) 



Hence we are left, in the approximation made, with a 
damped-oscillator equation of motion for the upper-state 
probability amplitude, where A(uj m ) and 5u> m , respec- 
tively, are given by Eqs. (58) and (M) [or Eqs. ( |60[ ) and 
(|6l|)1- Obviously, when uja = w m and 



ft 3> Su> r , 



(67) 



(i.e., strong coupling), then damped Rabi oscillations are 
observed: 



\C u {t)\ 2 = e- 5uJmt cos 2 (fit /2). 



(68) 



Note that in the opposite case where ft <C 5w m the solu- 
tion of Eq. ((65)) is \C u {t)\ 2 =e- At with A from Eq. @) 
for uja =w ra , which is in agreement with Eq. ([34|). Equa- 
tions (|66|) and (|67| ) together with Eqs. (|5^) and ( |59| ) [or 
Eqs. (|60|J and (|61|)1 provide us with an easy rule of thumb 
for deciding whether the strong-coupling regime is real- 
ized or not. 

In order to obtain the exact solution to the problem, 
we have also solved the basic integral equation (p8| ) [to- 
gether with the kernel function (B8j)] numerically. Typi- 
cal examples of the temporal evolution of the occupation 




FIG. 6. The temporal evolution of the occupation proba- 
bility |C u (i)| 2 of the upper atomic state is shown for c( = At, 
uja = 0.9999 ujt, A)At/(2c) = 10" 5 , and (a) lu p =3uj t , #2 = 
30.00197 At, (b) uj p = 1.5 wt, R 2 = 30.00179 A T . The solid 
Unes correspond to the exact solution and the dashed lines 
to the approximate analytical solution [solution of Eq. (|S5|)] 
shifted forwards by (a) A At =0.009 and (b) A At =0.02. 
For comparison, the exponential decay in free space is shown 
(dashed-dotted line). 



probability of the upper atomic state are shown in Fig. |5| 
for the case where the atomic transition is tuned to the 
resonance line in the center of the band gap. The figure 
reveals that with increasing value of the intrinsic damp- 
ing constant 7 of the wall material the Rabi oscillations 
become less pronounced, in agreement with Eqs. ( |5E} ) , 
(59), (|6fj), and the inequality ^6^). From Eqs. (58) and 
(59) it is seen that (in the approximation made there) the 
product A(uj m )5uj m ~c/ R2 does not vary with 7, whereas 
8u m linearly increases with 7. According to Eq. (|6^), 
f2ocy7 and thus fl/5uj m oc 1/^/7, i.e., the condition of 
strong coupling ( |67| ) becomes violated with increasing 7. 
Physically, increasing 7 means increasing probability of 
(irreversible) absorption of the emitted photon by the 
wall material and therefore reduced probability of (re- 
versible) atom-field energy exchange. It should be men- 
tioned that with increasing 7 the (very small) probability 
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that the photon (irreversibly) leaves the cavity also in- 
creases. Note that the variation of 7 in Fig. || leaves the 
imaginary part of the refractive index nearly unchanged, 
m ~ 1.2, while the (small) real part tir, slightly increases 
with 7 (see Fig. ||). 

The examples of the temporal evolution of the occu- 
pation probability of the upper atomic state shown in 
Fig. [| refer to the case where the atomic transition is 
tuned to a resonance line closest to lut below the band 
gap. According to Eqs. (|6(]), (|6l|), (|6rj), and the condi- 
tion (^37), the resonance frequencies close to wt are most 
favorable for realizing the strong-coupling regime in the 
range of normal dispersion, because of the rising real part 
of the refractive index [see Fig. ||(a)]. As expected, the 
strength of the the Rabi oscillations now varies with the 
plasma frequency wp such that they are less pronounced 
for small values of wp. Obviously, decreasing wp means 
increasing input-output coupling, i.e., increasing proba- 
bility that the emitted photon leaves the cavity instead 
of back-acting upon the atom. 

Finally, Fig. |^ presents a comparison between the exact 
solution, and the approximate analytical solution [solu- 
tion of Eq. (|65|)]. To compensate for the short-time inac- 
curacy of the analytical solution, it is somewhat shifted 
forward in time. The agreement between the exact solu- 
tion and the (shifted) analytical solution is quite good. 
Obviously, Eq. ( |65| ) predicts an initial decay somewhat 
faster than the exact one. In fact, the spontaneous decay 
is accelerated only gradually under the back-action of the 
radiated field being multiply reflected at the boundaries 
p3|p7t , and the single-resonance coupling is established 
only after a certain interval of time. The sharper is the 
cavity resonance, the shorter is this interval, which is 
fully confirmed by the figure. 



range within the band gap in a unified way. Whereas in 
the range of normal dispersion the cavity input-output 
coupling dominates the strength of the atom-field inter- 
action, the dominating effect within the band gap is the 
photon absorption by the wall material. 

In the study of the spherical-cavity problem, we have 
assumed that the atom is placed at the center of the 
cavity, which has drastically reduced the mathematical 
effort, because only the spherical Bessel function of or- 
der n — 1 contributes to the Green tensor. The price 
we paid is that the interaction of the atom with cavity 
excitations of large n, which correspond to high-Q whis- 
pering gallery modes and concentrate near the surface 
by repeated internal reflections has not been included in 
the analysis. Since the complete Green tensor is known, 
there is of course no obstacle to perform the calculations 
for an arbitrary position of the atom. In particular, when 
the atom is near the surface then nonradiative energy 
transfer from the atom to the absorbing medium can sub- 
stantially contribute to the process of spontaneous decay. 
Further investigations are also necessary in order to give 
a more detailed analysis of the evolution of the emitted 
radiation, to answer the question of the ratio of photon 
emission to nonradiative decay, and to extend the theory 
to multilevel atom-field interactions. 
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IV. SUMMARY AND CONCLUSIONS 

We have developed a formalism for studying sponta- 
neous decay of an excited atom in the presence of ar- 
bitrary dispersing and absorbing dielectric bodies. The 
formalism is based on a source-quantity representation 
of the electromagnetic field in terms of the Green tensor 
of the classical problem and appropriately chosen bosonic 
quantum fields. It replaces the standard concept of mode 
decomposition which fails for complex permittivity. All 
relevant information about the bodies such as form and 
intrinsic dispersion and absorption properties are con- 
tained in the Green tensor. It is worth noting that the 
Green tensor has been available for a large variety of con- 
figurations such as planarly, spherically, and cylindrically 
multilayered media ]|(J . 

We have applied the theory to the spontaneous decay 
of a two-level atom placed at the center of a three-layered 
spherical microcavity, modeling the wall by a Lorentz di- 
electric. The formalism has enabled us to study both the 
range of normal dispersion and the anomalous-dispersion 



APPENDIX A: THE HAMILTONIAN IN THE 
DIPOLE AND ROTATING WAVE 
APPROXIMATIONS 

The Hamiltonian (|l8|) can be rewritten as 

H = H F + H A + H AF , (Al) 
H F = Jd 3 rJ dwfejft(r,w)-f(r,w), (A2) 

^ = J2^ + IJ d 3 rp A (r)£ A (r), (A3) 

a 

Haf = - y —Pa • A(r a ) + / d 3 rp A (r)<p(r), (A4) 

where we have ignored the small A 2 term. Note that in 
the Coulomb gauge [p Q , A] = 0. For a neutral atom with 
the nucleus being positioned at r A , the atomic dipole 
operator reads 

= yq a {r a - r A ) = yq a r a , (A5) 
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and the first term in the interaction part of the Hamilto- 
nian -Haf takes the form 



^^p Q -A(r Q ) - - E 



if i 



Mr a) 



A(r A ), (A6) 



where the dipole approximation has been employed to 
replace A(r Q ) —> A(ta) and p Q = (m a / (ih))[r a , Ha] has 
been used. 

Now we restrict ourselves to a two-state model of an 
atom with upper state \u) and lower state |^). These are 
eigenstates of the unperturbed part of the Hamiltonian 
.Ha with the eigenvalues tiuj u and tkoi , respectively. Then 
one can write 



H A = huj u \u)(u\ + Hu>i\l){l\, 

\u)(u\ + \i)(i\ = i. 



(A7) 
(A8) 



In this atomic state space, the dipole operator fi A has the 
matrix elements (it |Aa 10 = (Z|/t A |w) = fi and M Aa \u) — 
(l\fi A \l) = 0. Using them and Eqs. @ - @, we 
arrive at 



V — p a ■ A(r Q ) 
£ — ' m, 



(a frt) 



du — — E (rA, oj) — H.c. 

UJ 



o-tE^(+)( rA ) .^t + H.c. 



(A9) 



where a = \l)(u\, = \u)(l\, lua — — ^i, and to = 
lua is set in the integral, because of the rotating wave 
approximation. 

In order to deal with the second term in Haf, we ex- 
pand /3a (r) in a multi-polar form and retain only the first 
non- vanishing term 

p A (r) ~ ^2q a 8(r - r A ) 

a 

-V • jff(r-rA)y)ga(r a - r A ) 



(A10) 



Then we have 

d 3 r p A (r)tp(r) 



-J rf 3 r{V-[ ( 5(r-r A )AA]}^(r) 

J d 3 r[8(r - r A )A A ] • V0(r) 
-Aa-E^pa) 



^eHW^a) /u + H.c. 



(All) 



where integration by parts and Eq. ( |15[ ) have been em- 
ployed for deriving the second and the third equation, 



respectively, and the rotating wave approximation has 
been used in de rivin g the forth equation. Combining 
(lAj), jA9|), and (|ATl|) gives 



^EW(r A ) /u + H.c. 



(A12) 



Equations flAip , (|A2|), (A12), and a subtraction of 
(h/2)(w u + u>i)I from H A , Eq. ( jA7| ), lead to the Hamil- 
tonian (^l|). 



APPENDIX B: THE GREEN TENSOR OF THE 
SPHERICAL CAVITY 



Following [|30|,f53| , we write the Green tensor of the cav- 
ity in Fig. |l| in the form 

G(r, r» = G v (r, r', u)8 fa + G^ (r, r», (Bl) 

where G v (r, r',a;) represents the contribution of the di- 
rect waves from the radiation sources in an unbounded 
medium, which is vacuum in our case, / and s denote 
the layers where the field point and source point locate, 
8f s is the usual Kronecker symbol, and the scattering 
Green tensor G^ 8 ' (r, r', uj) describes the contribution of 
the multiple reflection (/ = s) and transmission (/ =/= s) 
waves. In particular, G^ 13 \r, r', uj) and G ( - 33 ^(r, r', uj) 
read as t33|] 



e,o n— 1 m— 

J 2n + l (n-m)\ 
n(n+l) (n+m)\ 



+^ 3 (w)N g m (r, fci)N . nm (r' , fc 3 ) 

e.o n— 1 m=0 

J 2n + l (rc-m)! 
X \n(n+l)(n+m)! (2 d ° m) 

Cg(w)M. nn ,(r,*i)M. Bm (r' 1 *3) 



+C# (w)N - nm (r , fci )N . nm (r' , fc 3 ; 



(B2) 



(B3) 



where 



ki = k 3 



k 2 = \Muj) 



Ul 



and M and N represent TM- and TE- waves, respectively, 
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M >m (r, k)= T^Mkr)P™(cosd) ( Sm ) (m$e e 
sm (/ \ cos i 

■ n sdP.™(cos9) /cos\ . . 

N . nm (r, fc)= " (n . + 1) j» (fcr)F™ (cos 9) ( C ° S ) {m(p)e r 
kr V sm / 



1 d[rj„(fcr)] 



A' 7 



dP™(cos0)/cos\ 

(m0)e e 

\ sm/ 



M 



(B5) 



with j„(a;) being the spherical Bessel function of the 
first kind and P™(x) being the associated Legendre func- 
tion. The superscript (1) in Eq. (B2) indicates that in 
Eqs. (B4) and (B5) the spherical Bessel function j n (x) 
has to be replaced by the first-type spherical Hankel func- 
tion hn^ (x) 



The coefficients A^ N and Cf I N in Eqs. (jB2J) and 



arc defined by 



rpM ,N rpM ,N rpM ,N 

1 F1 1 F2 1 P1 



T, 



M,N 
PI 



rpM.N r/M,N r>M,N ' 
J-Fl tLpi -Kp2 



A 13 



T 



M,N 
P2 



R 



M,N 
P2 



R 



M,N 
PI 



T, 



M,N 
Fl 



T, 



M,N 
PI 



where 



R 



k f+ lH (f+i)f H ff 



k f H' ff H (f+1)f 



k f J' ff H {f+1)f 



R M 



k f+i J ff H {f+i)f 
k f+i J l f+i) f J ff ~ k f J 'ff J U+i)f 
Ff ~ k J+1 Jl 



R Pf - 



k f+1 H f 



f+1 n u+1)f n ff 



H' 



k f J (f+i)f H 'ff' 
k f H ff H ( f+ i)f 



k f+1 J' ff H {f+1)f - k f J ff H' u+1 



)f 



R N = k f+ l J (/ + l)/ J 'ff " k f J ff J (f+l)f 

Fi k f+i J (f+i)f H ff ~ k f J 'u+i)f H ff 



rpM _ 

x Pf - 



rpM 



rpN _ 
lp f - 



rnN _ 
1 Ff~ 

with 



fc /+i ( J (f+i)f H lf+i)f - J '(f+i)f H (f+i)f 

k f+i J ff H [f+i)f - k f J 'jf H U+i)f 
k f+i [j[f+i) f H(f + i)f - J(f+i)fH' (f+1 



)/ 



kf+iJ{ f+1)f H ff - k f J (f+1)f H' ff 
kf+i (j{ f+1 



i)f H (f+i)f 



J(f+i)fH' (f+1)f 



k f+1 J' ff H (f+1)f - kfJffH', 



'fl 

fc /+i ( J (f+i)f H [f + i)f 



(/ + !)/ 

J, (f+l)f H (f+vf 



k f+l J (f+l)f H 'ff - k f J if+l)f H ff 



Jil= jn(kiRl), 



(B6) 
(B7) 

(B8) 
(B9) 
(BIO) 
(Bll) 

(B12) 

(B13) 

(B14) 

(B15) 

(B16) 



H a =hW(k i R l ), 
p _ 1 d[pj n (p)} 
U ~ P dp 
, ld[ P h^{p)\ 



p=kiRi 



H 



P dp 



(B17) 
(B18) 

(B19) 



p=kiRi 



Note that A\ 3 N and C|| N are functions of n but not of 
to. When the atom is positioned at the cavity center, we 
have PI 



M-„ m (r,fc) 
N = nm (r,fc) 



fer->0 



At— 



-> (kr) 
(kr) n ~ 



(B20) 
(B21) 



In this c ase , only TM- waves with n — 1 contribute and 
the Eq. (B2) simplifies to 



G R EG< 33 '(r,r» 
Similarly, Eq. ( p32] ) reduces to 

G& 3 >(r,r>) 



ILO 

=r'—> 67TC 



(B22) 



= ^^(VlA^f.), (B23) 

isind d[rh^ (k^r)) 

r'— >0 Aitr dr 

(B24) 



0. 



(B25) 



APPENDIX C: DERIVATION OF EQ. Qsg) 

For an atom at the cavity center and z-oriented dipole, 
from Eqs. @ and (p23|) - (|B25[) we derive 



F r (r,r A ,L0 A ) = O(p- 2 ), 



(CI) 



F 9 (t,t a ,lu a ) = - 



k A n sin 9 
4ir 2 e p 



x / duj Im 
lo 



A^(Lj)e^/ c ] ((uj a -lu) + 0(p- 2 ), (C2) 



i^(r,r A , wa) = 0. 



(C3) 



Recalling the relatio n C( x ) — i/ ( x + *0), we perform the 
w-integration in Eq. (p2) to obtain 



duj Im 



{A^(w)e ,; ^/ c -c.c.} 



do; ■ 



2 Jq " cj - (cja + *0) 
-i 7 rA] v 3 (^A)e MA ' ,/c , (C4) 
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where we have (approximately) replaced A^(cu) by 
A]^(oja), extended the lower limit of the integral to — oo, 
and applied contour- i nteg ra tion techni que s. Combining 
Eqs. ©, ©, ©, @, @, and @, it is not dif- 
ficult to prove that 



T4 7 ~ flUJA 



\A$(ua)\ 2 
1 + ReCf (w A ) 



Taking into account that the free-space value Wq 
is obtained by setting A]$ 
yields Eq. (p3). 



1 and Cf : 



(C5) 



: htOA 



0, Eq. (05) just 
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